Wireless communications with diffuse waves 
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Diffuse, multiple-scattered waves can be very efficient for information transfer through disor- 
dered media, provided that antenna arrays are used for both transmission and reception of signals. 
Information capacity C of a communication channel between two identical linear arrays of n equally- 
spaced antennas, placed in a disordered medium with diffuse scattering, grows linearly with n and 
can attain considerable values, if antenna spacing a > A/2, where A is the wavelength. Decrease of 
a below A/2 makes the signals received by different antennas partially correlated, thus introducing 
redundancy and reducing capacity of the communication system. When the size of antenna arrays 
is well below A/2, the scaling of C with n becomes logarithmic and capacity is low. 
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Wireless communications in a disordered environment 
have recently attracted considerable attention jl], |2| [|. 
At first glance, scattering from randomly distributed 
heterogeneities disturbs the signal carried by the scat- 
tered wave (either acoustic or electromagnetic) and is 
expected to reduce the efficiency of communication. In- 
deed, the rate of error-free information transfer per Hertz 
of frequency bandwidth for a scalar, narrow-band com- 
munication channel between a transmitter and a re- 
ceiver is bounded by the channel information capacity 
= log 2 [1 + G\ /N] measured in bits per second per Hertz 
(bps/Hz) [|, [|. Here G — G(v) is the Fourier transform 
of the impulse response function G(t), describing the sig- 
nal at the receiver due to an infinitely short pulse emitted 
by the transmitter, and TV is the power of the Gaussian 
white noise at the receiver. Obviously, the decrease of 
the signal power |G| 2 due to scattering results in a re- 
duction of the channel capacity. However, in addition to 
the overall decrease of \G\ , scattering introduces random 
fluctuations of G in space. It was recently realized that 
one can make use of these fluctuations to overcome the 
reduction of C due to the drop of the signal power, pro- 
vided that the communication system contains a large 
number of antennas jl], This issue is very promis- 
ing for applications in wireless communications (mobile 
telephony in cities, indoor wireless local-area networks in 
buildings with complex structure, underwater communi- 
cations with acoustic waves, etc.). 

Assuming that the information about the scattering 
environment is available at the receivers (but not at the 
transmitters) and that uncorrelated noises at different 
receiving antennas have the same power N, the average 
capacity of a communication channel between arrays of n 
transmitting and n receiving antennas |L3[ can be defined 
as§ 

C = max(log 2 det[l„ + GQG+/A]), (I) 
Q 

where I n is the n x n unit matrix, G is a Green ma- 



trix (Giu gives the signal received by the receiver a 
due to the transmitter i), and Q is a non-negative def- 
inite covariance matrix describing correlations between 
the transmitted signals (with the constraint on the max- 
imum transmitted power TrQ < n). Angular brackets 
in Eq. ([!]) denote averaging over realizations of disorder. 
A rigorous analysis of the information capacity C in a 
disordered medium cannot be based uniquely on the ar- 
guments of the information theory ^, [s| and requires the 
physical understanding of scattering undergone by the 
waves carrying the information from transmitters to re- 
ceivers. In the present paper we analyze C in the frame- 
work of the model of diffuse multiple scattering, result- 
ing in a complicated, seemingly random spatial distribu- 
tion of scattered wave fields (so-called "speckles"). For 
the extreme cases of (a) a single speckle spot covering 
the whole antenna array (G has perfectly correlated en- 
tries) and (b) different antennas situated inside different 
speckle spots (entries of G are uncorrelated), we derive 
complete analytical expressions for C. We show that C 
increases during the continuous transition from the case 
(a) to the case (b) upon increasing antenna spacing (or, 
equivalently, upon decreasing correlations between the 
entries of G) and changes its asymptotic scaling with n 
from C oc Inn to C cx n. 

To be specific, we consider two identical linear arrays 
of equally-spaced transmitting/receiving antennas placed 
in a disordered medium. The distance L between the ar- 
rays is assumed to be much greater than the mean free 
path £ for waves in the medium, while the array size 
d = (n - l)a < I and d < {XL) 1 / 2 (Fresnel limit), 
where a is the spacing between adjacent antennas and 
A is the wavelength. In the majority of practically im- 
portant cases A <C t, and hence the propagation of waves 
from transmitters to receivers is diffusive . Entries G\ a 
of the matrix G can be then treated as complex Gaussian 
random variables with zero mean, equal variances, and 
possibly nontrivial correlations H 
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2 sin(fcArjj) sin(fcAr ai a) 
kArtj kAr a /3 



(2) 
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where k — 27r/A, Ar^ = \i — j\ a, and Ar a ^ = \a — j3\ a. 
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FIG. 1: Average information capacity per antenna of a com- 
munication channel between two arrays of n antennas in a 
disordered medium assuming statistical independence of the 
entries of the Green matrix G: results obtained using the ran- 
dom matrix theory (solid line) and Monte Carlo simulation 
(symbols) are shown for the signal to noise ratio S/N — 100. 
Horizontal dashed line shows the asymptotic value of C/n for 
n — * oo. Inset: Capacity per antenna as a function of S/N 
for single transmitting and single receiving antennas (dashed 
line) and an infinitely large number of antennas (solid line). 



The simplest cases to consider are (a) kd — > and 
(b) ka — run (m = 1,2,...). In the former case, 
all Gi a are perfectly correlated and one finds C = 

exp[l/(n 2 S/A0]£i[l/(n 2 SyA0]/ln2, wh ere E x is the ex- 
ponential integral function, S — no 2 is the average sig- 
nal power received at each receiver assuming independent 
signals from transmitters, and C oc Inn for n 2 S/N 3> 1. 
In the case (b), Q = I„ in Eq. (fi]) || and the machinery 
of the random matrix theory fllCfl can be employed for 
further analysis of the problem. In particular, it is useful 
to rewrite Eq. ([!]) as 



C=(£>g 2 [l + (S/iV)Mi 



(3) 



where \ii are the squares of the singular values of the ma- 
trix G+/S 1 / 2 with the joint probability density function 



JJ( Mi -^) 2 , (4) 



where A n is a normalization constant. The moment gen- 
erating function F[x) of a random variable c, defined 
by Eq. (|^) without averaging, is obtained by averaging 
exp(xc) with help of Eq. (0) O], and C = (c) is then 



C = T X F ^ 



x=0 



Tr (P^R) 



where P and R are n x n matrices: 



Pij = (i + i-2) 



(5) 



(6) 



du In [1 + [S/N)n] n i+ i~ 2 exp(-nu). (7) 



Equations provide an efhcient way for capacity 

calculation at arbitrary n and S/N (see the solid line in 
Fig. 0) and agree perfectly with the direct Monte Carlo 
simulation of capacity using Eq. (|) (symbols in Fig. |). 
Note that the growth of C with n is linear at large n and 
hence is much faster than in the case (a). 

An alternative way of calculating capacity consists in 
performing averaging directly in Eq. (|^) using Eq. (Q). 
This yields 



C/n 



d/ilog 2 [1 + (S/N)n]f n (n), 



(8) 



where = / °° d^ 2 ■ ■ ■ d^„ p n (^i, ■ ■ ■ , Mn) can be 

evaluated by a direct integration, at least at moderate 
n: /i(u) = exp(-^t), / 2 (/i) = 2exp(-2/i)[l - 2u + 2u 2 
/ 3 (m) = 3exp(-3u)[l - 6u + 18u 2 - 18/i 3 + (27/4)/z 4 
etc. The values of capacity obtained then from Eq. (^ 
coincide exactly with that following from Eq. (||). An 
asymptotic expression for f n (lJ>) at n > 1 can be found 
in the framework of the random matrix theory [l2| : 
fM = (2^)- 1 (4//i-l) 1 /2 for < n < 4 and /^(u) = 
otherwise. Eq. (H) then yields 

C/n = {S/N) 3 F 2 (1,1,3/2; 2,3; -4 5/iV)/ In 2, (9) 

where 3-F2 is the generalized hypergeometric function. 
This result is shown in the main plot of Fig. [l]by a dashed 
horizontal line. We find C/n ~ logjJ S/N) — log 2 e for 
S/N 3> 1 (see also the inset of Fig. [l]). It is worthwhile 
to note that C/n decreases monotonically with n, while 
the difference between the values of C/n at n — 1 and 
n > 1 never exceeds 7%. 

We now allow ka to take arbitrary values, thus in- 
troducing correlations between the entries of the Green 
matrix G. Eq. (Q) can again be reduced to Eq. (||) with 
/ii denoting the squares of the singular values of the ma- 
trix Q 1 / 2 G+/S' 1 / 2 , where Q is chosen to maximize the 
result. The joint probability density function of /x^, anal- 
ogous to Eq. (|j), is not known in this case, but one can 
still calculate C in the large-n limit [Q, |llj. The idea is 
to represent the moment generating function F(x) of a 
random variable c, defined by Eq. d3) without averaging, 
as a multiple Gaussian integral (the so-called "replica 
trick" ) and then to perform the integrations using saddle 
point methods in the limit n 3> 1. The maximization 
of C = (d/dx)F(x)\ x= o over Q is then accomplished by 
requiring SC < for all allowed small variations SQ of 
the optimal matrix Q. This yields a system of nonlinear 
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FIG. 2: Average information capacity per antenna of a com- 
munication channel between two identical linear arrays of n 
antennas in a disordered medium assuming antenna spacing 
a and S/N = 100. Solid lines correspond to n = 1 (exact 
result), n — 2 (Monte Carlo simulation) and n = 100 (asymp- 
totic large-n result). Dashed line shows the result obtained 
using the asymptotic large-n theory at n = 2. Inset: Capacity 
per antenna as a function of n for different antenna spacings a. 
Results for ka — n/2 and 7r/4 are asymptotic large-n results. 



equations for the eigenvalues of Q and some auxiliary 
variables that can be solved numerically. We refer the 
reader to Refs. [ll] and [l^ for an exhaustive account of the 
theoretical approach and the algorithm of the numerical 



calculation. The results that we obtained for identical 
linear arrays of equally-spaced antennas are presented in 
Fig. 0. As follows from the figure, at ka > ir correlations 
between the entries of G are too weak to affect C signifi- 
cantly and the latter remains very close to its maximum 
value, given by Eqs. (||) or (|j). In contrast, correlations 
become important when ka decreases below n, leading to 
a significant drop of C. The dashed line in Fig. || results 
from the asymptotic large-n theory |l], [ElJ |T^| with n = 2 
and is shown for illustration purposes only. Its closeness 
to the Monte Carlo result (solid line for n = 2) testifies 
a qualitative validity of the large-n theory even at small 
n. As ka increases, the scaling of C with n changes from 
C oc Inn at ka = to C oc n at ka > n (see the inset of 
Fig. H). We note that even at < ka < n, there is still 
a significant gain in capacity as compared to the case of 
ka = 0: e.g., at n = 100, C/n for ka — n/2 (tt/A) is 
almost 20 (13) times larger than for ka = 0. 

In conclusion, we have studied the information capac- 
ity C of a wireless communication channel in a disordered 
medium, assuming multiple diffuse scattering of waves 
that carry information, and taking into account fluctua- 
tions of wave fields in space (speckles). Although multi- 
ple scattering reduces the received signal, it allows for a 
dramatic increase of capacity in the case of communica- 
tion between two antenna arrays, provided that antenna 
spacing a > A/ 2, where A is the wavelength. Namely, 
for identical linear arrays of n equally spaced receiv- 
ing/transmitting antennas, scaling of C with n changes 
from C oc Inn for kd <C tt to C oc n for ka > 7T, where 
d = a(n — 1) is the array size and k = 2n/ A. Even at 
< ka < 7r an important gain in capacity is possible as 
compared to ka = 0. 

The author is indebted to Prof. R. Maynard for helpful 
discussions and critical reading of the manuscript. 
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For simplicity, we consider the transmitting and receiving 
arrays to consist of the same number n of antennas. 



